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Abstract 

A relation bet-ween semi-direct sums of Lie algebras and integrable couplings of 
lattice equations is established, and a practicable -way to construct integrable cou- 
plings is further proposed. An application of the resulting general theory to the 
generalized Toda spectral problem yields f-wo classes of integrable couplings for the 
generalized Toda hierarchy of lattice equations. The construction of integrable cou- 
plings using semi-direct sums of Lie algebras provides a good source of information 
on complete classification of integrable lattice equations. 
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1 Introduction 

Integrable couplings have been receiving gro-wing attention recently. A few -ways to con- 
struct integrable couplings are presented by using perturbations [HEIIS], enlarging spectral 
problems jHE], and creating new loop Lie algebras [EllZI- 

The problem of integrable couplings can be expressed as follows 0: For a given in- 
tegrable system, how can we construct a non-trivial system of differential equations which 
is still integrable and includes the original integrable system as a sub-system? Obviously, a 
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change of orders of equations in a system does not lose integrability of the system. There- 
fore, up to a permutation, an integrable couphng of a given integrable system Ut = K{u) 
is given by a bigger and triangular system: 

Ut = K{u), vt = S{u,v). 

The vector- valued function 5* should satisfy the non-triviality condition dS/d[u] ^ 0, where 
[u] = {u, D^u, D^u, ■ ■ ■ ) and D^u denotes a vector consisting of all derivatives of u of order 
n with respect to the space variable x. The above non-triviality condition means that the 
second sub-system involves the dependent variables of the first sub-system (i.e. the original 
system), and thus it guarantees that trivial diagonal systems with S{u,v) = S{v) are not 
within our discussion. 

A basic integrable coupling of an integrable system ut = K{u) is given by 

Ut = K{u), Vt = K\u)[vl (1.1) 

which can be generated by a perturbation around a solution of the system ut = K{u) 1 . 
In the above system and elsewhere throughout this paper, P'{u)[v\ denotes the Gateaux 
derivative of P{u) = P{u, D^u, ■ ■ ■ ) with respect to m in a direction v, i.e., 



d 



d 



P'{u) [v] = —P{u + ev) = —P{u + ev, D^u + sD^v, ■ ■ 



e=0 



e=o de 

Obviously, the second sub-system Vt = K'{u)[v] in the above integrable coupling ()1.1|) is 
linear with respect to v. Moreover, a symmetry S{u) of the system ut = K{u) leads to 
a solution {u,S{u)) to the integrable coupling (jl.lj) . However, the second component v 
of a solution (-u, v) to the integrable coupling (jl.lj) is generally not a symmetry of the 
system ut = K{u). This is because v satisfies the linearized system vt = K'{u)[v] only 
for one solution, not for all solutions of the system Ut = K{u). Therefore, the simple 
integrable coupling (jl.lj) is already a generalization of the symmetry problem. Another 
basic integrable coupling of an integrable system Ut = K{u) reads as 

Ut = K{u), Vt = K'{u)[v] + K{u). (1.2) 

This system has a set of hereditary recursion operators [2j 

(3Mu) 

/3i<I>'(m)H+/52$(u) /3i$(n) 

with two arbitrary constants /3i and (32-, if the original system ut = K{u) has a hereditary 
recursion operator $(?i). Therefore, integrable couplings possess richer integrable structures 
than the original integrable systems. 
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The study of integrable couphngs provides clues towards complete classification of 
integrable systems. Let us first observe classification of square matrix spectral problems 
through the Jordan blocks under similar transformations of matrices. Each triangular Jor- 
dan block corresponds to an undecomposable sub-system in a given integrable system. Now, 
note that an arbitrary Lie algebra has a semi-direct sum structure of a solvable Lie algebra 
and a semi-simple Lie algebra [H], and we will see that semi-direct sums of Lie algebras can 
result in integrable couplings. These imply that the study of integrable couplings through 
semi-direct sums of Lie algebras is an inevitable step towards complete classification of 
integrable systems with an arbitrary number of components, from a point of view of Lie 
algebras. 

The study of integrable couplings also generates interesting mathematical structures 
such as Lax pairs with several spectral parameters jHlEUE], integrable constrained fiows 
with higher multiplicity fTH, I12j. local bi-Hamiltonian structures in higher dimensions ^3] 
and hereditary recursion operators of higher order [21 E] ■ Very recently, we have proposed 
a relation between semi-direct sums of Lie algebras and integrable couplings of continuous 
soliton equations, which provides an interesting and systematic approach to integrable 
couplings of continuous soliton equations [T3]. In this paper, we would like to discuss the 
problem of discrete integrable couplings and develop a theory for constructing discrete 
integrable couplings by use of semi-direct sums of Lie algebras. 

Throughout our discussion, we denote by E the shift operator, write 

{E"'x){n) = x^"''\n) = x{m + n), where x : Z ^ R, m,n G Z, (1.3) 
and adopt an inverse of the difference operator ii^ — 1 as follows 

—1 oo 

iE-ir = l{Y: -E)^'- (1-4) 

fc=— oo k=0 

Let G be a matrix Lie algebra with the standard Lie bracket [A, B] = AB — BA, and closed 
under matrix multiplication: AB G G for all A,B E G. We assume that an integrable lattice 
equation (or system) of evolution type 

ut = K{u) = K{u,Eu,E-\---) (1.5) 

is associated with G, where u = u{n,t) is a dependent variable. More precisely, there is a 
pair of square matrices U and V in G, called a Lax pair, so that the discrete spatial matrix 
spectral problem 

E(f) = U(f) = U{u,X)(j) (1.6) 
and the associated discrete temporal matrix spectral problem 

(j)t = V(f) = V{u,Eu,E-^u,--- ■,X)(f), (1.7) 
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where A is a spectral parameter and is an eigenfunction, generate [T6t ITTj the integrable 
lattice equation p.5|) through their isospectral (i.e., Xt = 0) compatibility condition 



Ut = {EV)U - UV, 



(1.8) 



which is called a discrete zero curvature equation. In other words, we have 



U'{u)[K] = {EV)U-UV, 



(1.9) 



where U'{u)[K] denotes the Gateaux derivative as above. In a non-isospectral case, e.g., 
Xt = /(A), then we have 



where U\ is a partial derivative of U with respect to A. Based on (jl.fij) and p.7|l . the lattice 
equation ()1.5|1 can often be solved by the inverse scattering transform (for example, see J18j). 
There are also a few interesting Lie algebraic structures hidden behind the equation ()1.10j) 
(see jTTj for more information). An integrable hierarchy and its master symmetry hierarchy 
usually correspond to the isospectral case and the non-isospectral case Aj = A", n > 0, 
respectively. These two hierarchies constitute a semi-direct sum of Lie algebras, each of 
which consists of symmetries in one hierarchy. The spatial matrix spectral problem (jl.fij) 
is our starting point in constructing discrete integrable couplings. The closure property of 
the Lie algebra G under matrix multiplication guarantees that {EV)U — UV is still in G 
so that the discrete zero curvature equation ()1.8|1 makes sense. 

In what follows, we are going to establish a relation between semi-direct sums of Lie 
algebras and integrable couplings of lattice equations and a technically-practicable way 
to generate integrable couplings through semi-direct sums of Lie algebras. The resulting 
general theory will be used to generate two classes of integrable couplings for the generalized 
Toda hierarchy presented in JHl- It will also be indicated that the construction of integrable 
couplings using semi-direct sums of Lie algebras provides a good source of information 
about classification of integrable lattice equations. A few concluding remarks will be given 
in the last section. 

2 Constructing integrable couplings using semi-direct 
sums of Lie algebras 

2.1 Generating scheme 

Assume that the lattice equation ()1.5|) has a Lax pair {U, V) in a matrix Lie algebra G 
closed under matrix multiplication. 
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To construct an integrable coupling of the lattice equation ()1.5|1 . we use semi-direct 
sums of Lie algebras to enlarge the original Lie algebra G. Take another matrix Lie algebra 
Gc closed under matrix multiplication and then form a semi-direct sum G of G and Gc- 

G = G € Ge. (2.1) 

The notion of semi-direct sums means that G and Gc satisfy 

[G,Gc]CGc, (2.2) 

where [G, GJ = {[A^B]\A E G, Be Gc}. Obviously, Gc is an ideal Lie sub-algebra of 
G. The subscript c here indicates a contribution to the construction of couplings. We also 
require that the closure property between G and Gc under matrix multiplication: 

GGc, GcG 'i^ Gc, (2-3) 

where G1G2 = {AB \ A G Gi, B G G2}, to guarantee that a Lax pair from the semi-direct 
sum G can generate a coupling system. Note that the two different binary operations were 
used in the above closure properties in ()2.2|) and ()2.3|) . 

Now choose a pair of new Lax matrices in the semi-direct sum G of Lie algebras: 

U = U + Uc, V = V + Vc, t/c, K e Gc, (2.4) 

and make a pair of enlarged discrete spatial matrix spectral problems 

E^ = U^ = U(u,XU, 

_ _ _ (2.5) 

4>t = V(p = V{u,Eu,E~^u--- ■,X)c^, 

where the matrix Uc in U introduces additional dependent variables and u consists of both 
the original dependent variables and the additional dependent variables. In addition, the 
matrix Uc could depend on the spectral parameter A, and the matrix in really does 
almost in all cases. Based on the closure properties of G, Gc and between G and Gc, it is 
easy to see that 

{EV)U -tJV = [{EV)U - UV] 

+{[iEV)Uc - UcV] + [{EVc)U - UVc] + [{EVc)Uc - f/cK]} e G E Gc- 
Therefore, under Ut = K{u), the corresponding enlarged discrete zero curvature equation 

Ut = {EV)U - UV (2.6) 

precisely presents 

Ut = (EV)U - UV, 

(2.7) 

Uc,t = [{EV)Uc - UcV] + [{EVc)U - UVc] + [{EVc)Uc - UcV]. 
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The first equation above is equivalent to the lattice equation ()1.5|1 . and hence, this is a 
coupling system for the lattice equation ()1.5|) . 

The whole construction process above shows that semi-direct sums of a given Lie 
algebra G with new Lie algebras provide a great choice of candidates of integrable couplings 
for the lattice equation ()1.5|) generated from the Lie algebra G. 



2.2 Realizations by particular semi-direct sums 

To shed light on the above general scheme of constructing coupling systems, let us introduce 
the following particular class of semi-direct sums of Lie algebras: 



G = G e a, G 



r /o 



diag( A 




Gr 



\ 1 



5,; 



B 



At+l,At+l 



I v 



B 



(2.8) 



J ) 



where A, Ba are arbitrary square matrices, A is of the same order as U and the partitions of 
matrices in G and Gc are the same. Obviously, Bij, j < fi, are square but Bij, j > fi+l, may 
not; and all closure conditions of G, Gc and between G and Gc under matrix multiplication 
are satisfied. 

Define 



Ug := diag(2^_;^, 0^_^), Vg := diag(y^_^, 0^_^). 

Note that Ug and Vg in G generate the same lattice equation as U and V, and thus for 
integrable couplings, the corresponding enlarged spectral matrices U and V in the semi- 
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direct sum G ^ Gn can be chosen as 



U = Ug + Ug,c 



\ 



u 



u, 



V = Vg + Vg,c 



\ 



V 



K+1 



\ 



where the first two matrices Ug and Vg play the {U, y)-part and the second two matrices 
Ug,c and Vg,c play the (f/c, Vc)-part in the pair of U and V defined in ()2.4|1 . It is not difficult 
to see that the resulting coupling system ()2.7|) becomes 



' Ut = {EV)U - UV, 

= (-£'Vi+i,i+i)f/j+i^i+i — f/j+i^i+i^+i^j+i, fi < i < ly, 

j 



(2.9) 



k=i 



where Uu = U and Vu = V, I < i < fi. 
In particular, first, if we take 



U 



V 







fV Vi2 ■■■ Vi,^^i\ 



V 



V V, 



V 







then the coupling system ()2.9|) becomes 

Ut = {EV)U - UV, 
j 

Uij,t = J2i(^^^k)Ukj - UikVkj], I <i< 3 < fi + l, 



(2.10) 



k=i 
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where Uu = U, Vu = V, I < i < jj,, and f/, 

u ■•• I 



u 



Vii = 0, i = jj, + 1. Second, if we take 

fV ■■■ 
V ■•• I 



V 



\o uj 

then the couphng system ()2.9|1 becomes 

f/t = {EV)U - t/V, 



\0 



ai 



5^ [{EV^,)U^,-UM, i<^< 



fc+i=j,fc,i>0 



where f/„„ = U and K,, 



(2.11) 



We remark that here we have just presented one class of semi-direct sums of Lie al- 
gebras, together with two specific examples. It is interesting to construct other possible 
realizations, especially those which could carry essential information for keeping integrable 
properties of the original lattice equations. 



2.3 Linearly dependent case on the spectral parameter 

Let us now assume that the spatial spectral matrix U depends linearly on the spectral 
parameter A (see, for example, [El UHl EDI EI] ) : 

U = Uiu, A) = XUo + U,, ^ = ^ = 0. (2.12) 

Consider two specific examples of the enlarged spatial spectral matrices introduced in the 
last sub-section: 

'u f/„\ fu u,\ au„ 



Note that the sub-matrices Ua in the above two enlarged spatial spectral matrices could 
be of different sizes. As in the continuous cases IHHS], suppose that 



^0 / V w 

with 

W = Y: W.X-\ W^=Y: ^ = 0' ^ = 0' (2-14) 

j>0 j>— no 

where no > is a proper integer, solve the corresponding enlarged discrete stationary zero 
curvature equations 

{EW,)U,-UM = 0, 2 = 1,2, (2.15) 
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respectively. 

Then for each m > 0, choose 



't/h t/["1\ /a a 



y V 

where A^ and Am,a do not depend on A and satisfy 

(EA^)f/o - t/oA„ = 0, f/oA„,„ = 0; (2.16) 

and choose 



T/m T/^'^U /A A 



m 



where A^, and Am,a do not depend on A and satisfy 

(^A„)f/o - f/oA^ = 0, (^A^,a)f/o - UQ\n,a = 0. (2.17) 

The subscript + above denotes to select the polynomial part in A. Based on ()2.1()|1 and 
()2.1H) and using ()2.15j) . we can directly show that the enlarged discrete zero curvature 
equations 



namely, 



present 



and 



\ M and Ua,t^ = (EVH)f/„ + ^EV^^)u 



Ut^ = (A„), + [Uo, W^+i] - [f/i, A„], 



(2.18) 



Ut^ = (A„), + [Uo, Wm+i] - [Ui, Am], 

Ua,t^ = UoWa,m+l - iEWa,m+l)Uo (2.19) 
+ {EAm)Ua - UaArn + {EArn,a)Ul - UiAm,a, 

respectively. 

We remark that these two enlarged hierarchies in ()2.18|) and ()2.19|) share the enlarged 
discrete spectral problems 

E(j) = Ui(j), E(j) = U2(p, 
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respectively. Thus, all lattice equations in each of the two enlarged hierarchies can possess 
infinitely many common conserved densities except the original ones (see ^Hl |23 123 for 
a few concrete examples). Moreover, one can construct a specific non-degenerate bilinear 
form on G with the invariance property, to present Hamiltonian structures of the enlarged 
lattice equations by a generalized trace identity. The detailed analysis on those integrable 
properties will be left to a future publication. 

To sum up, each system of lattice equations in the hierarchy ()2.18|) or (j2.19|l can 
provide an integrable coupling for its first sub-system of lattice equations. In the next 
section, we will only discuss two examples of constructing enlarged lattice hierarchies, in 
the generalized Toda case presented in |19j. 



3 Integrable couplings of the generalized Toda hier- 
archy 

3.1 The generalized Toda equations 

Let us here recall the generalized Toda hierarchy ^H]- The corresponding discrete spatial 
spectral problem reads 

E(P = U{u,X)<p, U{u,\)={ ° ^ I, (3.1) 

where A is a spectral parameter, and a and (3 are two arbitrary constants satisfying a^+(3'^ ^ 
0. When a = and [3 = —1, ()3.1|) becomes the Toda spectral problem [T^ . 

Its stationary discrete zero curvature equation 

{EW)U -UW = (3.2) 



has the solution 



with 




a = ^aiX \ b = ^biX \ c=^CiA \ 

i>0 «>0 i>0 

where the coefficients are defined by the initial conditions: 

ao = ho = 0, Co = 0, 
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and the recursion relation: 

' Q+i - rhf^^ = 0, 



feffi + sb^^ + (a^^ + a.) = 0, 
(a-+\ - Oi+i) + s{a^p - ai) 
[ +a(r6,+i - cffi) + f3{rb, - cf'>) = 0, 



i > 0, 



(3.4) 



which are all difference polynomials in u with respect to the lattice variable n. Under the 
initial-value conditions 

ai|tj=o = Ci|„=o = 0, ai\u=o = bi\u=o = Cj|„=o = 0, i >2, 

the recursion relation ()3.4j) uniquely determine the lattice functions aj,6j and c,, i > 1. 
The first few lattice functions are 

ai = ar, bi = 1, Ci = r; 

a2 = —a'^r^-^h' — a^r^ — a'^rr^'^^ — ars — ars*^""^-* + /5r, 
62 = —a?" — cer^"^^ — s^^"'^-', C2 = — rs — ar^ — arr*^-*^). 



As usual, choose that 

Then it follows from ()3.4|1 that 

{EVjU-UVm 
Take a modification 



(A™a)+ (A™6)4 
(aA + /5)(A™c)+ -(A™a) 



m > 0. 







(1) 



A. 



brn+l 





and define the temporal spectral matrices 

=Vm + Am, m> 0. 







Then, a direct calculation leads to the following matrix: 



(aA + f3){cm+i - rhm+i) f3{cm - rbm) - s(a.l^^ - a„ 
This is consistent with [/f^, and thus, making the evolution laws 



(3.5) 



(3.6) 



(3.7) 
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the compatibility conditions 

Ut^ = {EV^^"^)U - UV^'^\ m > 0, 

of the discrete spatial spectral problem and the associated discrete temporal spectral 
problems ()H.7|1 give rise to the following hierarchy of lattice equations 

m>0. (3.8) 

This generalized Toda hierarchy is Liouville integrable ^Hj, and its Hamiltonian structure 
leads to infinitely many conservation laws and symmetries for every system in the hierarchy. 

Obviously, the first nonlinear lattice equation in the hierarchy is 

= r(s^^^^ — s) + ar(r^^^^ — r^^^), 

(3.9) 

= as{r — r^-^') + (3{r^'^' — r). 



When a = and P = —1, ()3.9p becomes the Toda lattice equation 

= r(s("^) - s), St,=r- r^; (3.10) 
and when a = 1 and /3 = 0, ()3.9|1 becomes the following lattice equation presented in |26j : 

= r(s("^) -s)+ r(r("^) - r^^^), St, = s(r - r^). (3.11) 

The lattice equation 1)3.111) is linearly independent of the Toda lattice equation ()3.1()|1 . 
There exist a voluminous literature on the Toda lattice equation, and its generalizations 
and solution structures (for example, see [211- |S2])- 

3.2 Integrable couplings from specific semi-direct sums 

The generalized Toda spectral problem ()3.ip linearly depends on the spectral parameter 
A, and thus we can write 

U=\ ° ^ ] =UoX + Ui, Uo= \ ^ °l,f/i=|° ). (3.12) 

We will also see that there is a difference between the two cases of a = and a 7^ in 
computing integrable couplings. 

Let us first consider the semi-direct sum of Lie algebras of 3 x 3 matrices: 



'A 0^ 




A G C[A,A-^] 0M2X2 L G, 



'0 





5gC[A,A-1]®M2xi 
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where C[A, A ^] Mmxn = span{ X'^A \ k E A E Mmxn }• In this case, Gc is an Abehan 
ideal of G € Gc- We define the corresponding enlarged spatial spectral matrix as 

U = Uiu, A) eGe Gc, Ua = Ua{v) = (^^^ , (3.13) 

where vi and V2 are new dependent variables and 

V = {vi, V2f, u = {u^, v^y = (r, S, Vi, V2f. 

Upon setting 

where is a solution to {EW)U — UW = 0, defined by ()3.3p . the corresponding enlarged 
discrete stationary zero curvature equation {EW)lf — UW = becomes 

{EW)Ua-UWa = 0, (3.14) 

which is equivalent to 
namely, 

{aX + (3)re={aX + (3)c^^hi-a^^^V2-{X + s)f. 
This system determines a solution for e and / as follows 

e= ^ e,X-\ f = J2fiX\ 

i>— no i>0 

where no = 1 if a = and no = if a 7^ (see ()2.14|1 for introduction of no). Now define 
the enlarged temporal spectral matrix as 

^''"l = , = {X^Wa)+ + A„,„ m > 0, (3.16) 

V J 

where V^^'^ is defined as in ()3.(ij) . To satisfy ()2.1f)|l . choose Am,a as 

(hm \ 
] , hm - arbitrary, m > 0. (3.17) 
-arhmj 

Then based on ()3.15p . we can compute that 

^-acj^l^it;! + are„+i + Z^+i/ V / \(3rhm - arsh 

-ac^^^^vi + arcm+i + fm+i - (3rhm + arshm^ 

13 



m > 0. 



Therefore, the m-th enlarged discrete zero curvature equation 
leads to 



-ac^^\^vi + arcm+i + fm+i - f3rh^ + arsK 



(3.18) 



together with the m-th generalized Toda equation in ()H.8|) . Therefore, we obtain a hierarchy 
of coupling systems defined by ()2.18|1 : 



Ut 



u 

V 



Krniu) 

Kra{u) = I , I , m > 

bm[u,v) 



(3.19) 



for the generalized Toda hierarchy ()3.8|) . 

Let us second consider the semi-direct sum of Lie algebras of 4 x 4 matrices: 



'A 
A 



A e C[A,A'^] 0M2X2 } , Gc 



'0 





2x2 



In this case, Gc is an Abelian ideal of G € Gc, too. We define the corresponding enlarged 
spatial spectral matrix as 



U = U{u,\) 



U f/„ 



eG^Gc, Ua = Ua{v) 



Vl V2 



V3 V4, 



t/ 

where Vi, 1 < i < 4, are new dependent variables and 



(3.20) 



If we set 



'W Wa\ e f 

W=\ , Wa = Wa{u,X) = 

W \g -e^ 



where W is a. solution to {EW)U — UW = 0, defined by ()3.3|) . then the corresponding 
enlarged discrete stationary zero curvature equation {EW)U — UW = becomes 



{EW)Ua + {EWa)U - UWa - UaW = 0, 



(3.21) 



which is equivalent to 

(e(i) + e) + (A + s)/(^) + (a(i) + a)v2 + b^^^v^ - bv^ = 0, 

-(aA + /5)r(e(i) + e) - (A + s)g + (aX + P)ic^^hi - cv^) - (a^^) + a)v3 = 0, (3.22) 
^(1) - (A + s)(e(i) - e) - (aA + f3){rf - c^^hi) - (a^^) - a)v4 - bv^ = 0. 
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This system can determine a solution for e, / and g as follows 

i>0 i>0 i>0 

Now, we define the enlarged temporal spectral matrix as 

\ i^ny 

where V^"^'^ is defined as in ()3.6|) . To satisfy p.l7|l . choose Am,a as 

\ , hm - arbitrary, m > 0. 
-arhra y 

Then based on ()3.22p . we can compute that 
(EV^H)f/^ + {EVl"'^)U - f/yj"' - f/a^'""' 

(1) Ml) I ^ (1) f I (1) 

-bm+iV3 J \r{as-(3)hm -arWhg\ 

,acj^^it;i - ac„,+iV4 - ar{e^^\^ + e„+i) - 5(^+1 - fe^+i^^s + '^(as - (3)hr, 

Jm+l + ''m+l^S + rim \ 

"^+1^2 - - {e^^\^ - e„,+i) - ar^^'^h^^ J 

Then the m-th enlarged discrete zero curvature equation 

leads to 

/ -aCm+iV2 + arfi^l^ + (fej^^i - 6^+1 + a^^/im \ 
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together with the m-th generahzed Toda equation in ()3.8|1 . Therefore, we obtain a hierarchy 
of couphng systems defined by ()2.19|) : 

ut^=( ] = KM =[ rr. r 1 , m > (3.26) 
for the generahzed Toda hierarchy ()3.8|1 . 

3.3 Illustrative examples 

We now work out two concrete examples as follows, one in each of the two above cases. 

Case of a = 1 and /3 = 0: Let us first compute an example of the hierarchy ()3.19|) . 
Assume that a = 1 and /? = for convenience, which corresponds to the lattice hierarchy 
presented in jSEI- In this case, we have no = in ()2.14j) . 

It directly follows from ()3.15|1 that 

(1) (1) f f 

where z > and reo = Cq^"^! ~ fo- We can then obtain that 
/o = -^Vu fi = r'^^^vi + V2, 

/2 = -[r(2)r(i) + (r(i))2 + r^^V + r^^^s^^^ + r^^'>s]vi - (r^^) + r + s)v2; 
reo = 1^1, rei = ^svi — |f2, re2 = r^^Vf i + rv2. 

If we choose 

h = 2^/o/i = -^viir^^'^vi + V2), ^ = const., 
then the vector-valued function Si defined by ()3.18|) becomes 

+ r + s)vi — ^rvi{A^'^Vi + V2) 
—r^^^svi — {r^^^ + s)v2 — ^rsvi{r^^^Vi + V2) 



Si{u,v) 



Therefore, the integrable coupling of the generalized Toda lattice equation ()3.1H) . defined 
by ()3.19|) . reads as 

rt-^ = r(s^^^^ — s) + r(r^^^^ — r^^^), 
St-, = s(r — r^^^), 

(3.27) 

Vi^ti = — (r^-^) + r + s)vi — C^rvi{r^^^Vi + V2), 
. "^2,41 = —r^^^svi — (r^^^ + s)v2 — ^rsvi{r^^^Vi + V2), 
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the second sub-system of which is nonhnear with respect to both sub-sets of dependent 
variables when ^ 7^ 0. 

Case of a = and [3 = 1: Let us second compute an example of the hierarchy ()3.26|) . 
Assume that a = and (3=1 for convenience, which corresponds to the Toda lattice 
hierarchy. In this case, we have no = 1 in ()2.14|) . 

We take the initial set of functions as follows: 



fo = 9o = 0, Co 
Obviously from ()3.22|) . we can have 



Ji+1 


= -sfi 


- {ef^ + e.) 


- (af ^ 


9i+i 


= -SQi - 


r(ef + Ci) 


- {4' 




- ei+1 = 




+ 4' 



where i > 0. It then follows that 



r /i = i 



,(-1) 
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r + V3, ei = V3 + r(l + ^^); 



/2 = + r)(l + vi''^) - - v^'^ - r(-i)(l + v!f^^^ -^"'^ 



- V) 



(-1) 



92 = 



= —s{r + vs) — r[t>3 -|- f 3 -|- r(-^)(l -|- V2) + r(l -|- t>2 )] + r^'^^Vi — rv^ 
— 62 = (1 — s){r^^^ — r + fg""""* — f 3 + r'^'^^V2 — ^^)- 



We can use the inverse formula p.4|) to compute 62 here, but as we will see, this is not 
necessary for computing the corresponding integrable coupling. 

Now if we choose 

hi = TjViei = ?7fi[f3 -|- r(l -|- t>2 "'^^)], rj = const., 
then the vector-valued function Ti defined by ()3.25p becomes 

-SV2 - (s + r(i))(l + V2) - r(l + v'f^^) + 



(1) 



s)w3 rs r[f5^^ + ^^3 + r(^)(l + V2) + r(l v!f^^)] 
_l_ 7-^^ _ rjrvilvs + r(l -|- "^2^"'"^)] 



V 



l)(r 



(1) 



,(1) 



r -|- f 3"' — f 3 -|- r^^)f2 — rv. 
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Therefore, the integrable couphng of the Toda lattice equation ()3.in|l . defined by ()3.26|1 . 

reads as 

r^j = r(s^^^^ — s), St-i = r^^^ — r, 

V2fy = -SV2 - (s + r(^))(l + V2) - r(l + v!f^^) + f 1 - fg^^ 
< -vs - V4 + 7]v'i^\v^^^ + r^^^l + V2)], (3.28) 

vsfy = (s^^^) + s)vs + rs + r[v^^^^ + V3 + r'-^^l + V2) + r(l + f^"^^)] 
—j-Wy^ + rv4 — rjrvilv^ + r(l + "'"^)], 

■^4,*! = (s — l){r^^^ — r + v^^^ — "^3 + r^^'^V2 — rt'g""'"'*), 
the second sub-system of which is nonhnear with respect to both sub-sets of dependent 
variables when 77 7^ 0. 

4 Conclusions and remarks 

A feasible approach to construct integrable couplings of discrete soliton equations has 
been proposed by taking advantage of semi-direct sums of Lie algebras, and the resulting 
theory has been applied to the generalized Toda hierarchy of lattice equations to generate 
integrable couplings for the hierarchy. The key point in our generating scheme is to establish 
a relation between semi-direct sums of Lie algebras and integrable couplings of discrete 
soliton equations. The underlying discrete matrix spectral problems are generated from 
semi-direct sums of Lie algebras, and the discrete Lax spectral matrices associated with 
given soliton equations play the non-ideal part in the semi-direct sums. 

In our analysis of the two specific semi-direct sums, we have seen that there is always 
an arbitrary modified term An,a- This indicates that higher order matrix spectral problems 
have more degrees of freedom in generating integrable systems. On the other hand, in all 
additional spectral sub-matrices such as Ua^ , one can take their dependence on the spectral 
parameter into consideration, and this will bring much more diverse integrable couplings. 

There are also other questions about integrable properties of the resulting enlarged 
lattice equations, even in the case where additional spectral sub-matrices are independent 
of the spectral parameter. For example, can we solve the enlarged lattice equations by the 
inverse scattering transform? The class of Lie algebras in ()2.8|) provides a few realizations of 
semi-direct sums of Lie algebras. Other possible realizations are still interesting, especially 
those which could carry significant information about integrable properties. Reductions of 
the presented cases of semi-direct sums, which keep the uniqueness property of discrete 
spectral problems (see [H]), could be good examples. 
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We would especially like to emphasize that we have been considering the problem of 
integrable couplings and the key is semi-direct sums of Lie algebras. The initial Lie algebras 
G associated with given integrable systems in our construction can be simple (e.g., see 
|29t l3Uj). but semi-direct sums of Lie algebras G are normally non-simple (see ^3]). Our 
examples in Section 2 are all non-simple, since the Killing forms on those semi-direct sums 
of Lie algebras G are degenerate. However, there still exist specific non-degenerate bilinear 
forms on the Lie algebras G, with nice invariance properties, and their corresponding 
generalized trace identities, which present Hamiltonian structures of the enlarged lattice 
equations. 

To conclude, semi-direct sums of Lie algebras provide a good source of matrix spectral 
problems for generating integrable systems, and thus the study of integrable couplings 
using semi-direct sums of Lie algebras will enhance our understanding of classification of 
integrable systems. We are expecting to see more research on related topics. 
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